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Abstract. We study both theoretically and experimentally the free induction decay 
(FID) of the electron spin associated with a single nitrogen-vacancy defect in high 
purity diamond, where the main source of decoherence is the hyperfme interaction 
with a bath of 13 C nuclear spins. In particular, we report a systematic study of the 
FID signal as a function of the strength of a magnetic held oriented along the symmetry 
axis of the defect. On average, an increment of the coherence time by a factor of a/5/2 
is observed at high magnetic field in diamond samples with a natural abundance of 13 C 
nuclear spins, in agreement with numerical simulations and theoretical studies. Further 
theoretical analysis shows that this enhancement is independent of the concentration 
of nuclear spin impurities. By dividing the nuclear spin bath into shells and cones, we 
theoretically identify which nuclear spins are responsible for the observed dynamics. 
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1. Introduction 

Color centers in solids have emerged as good candidates for quantum information 
processing as they provide optical access on demand to a quantum degree of freedom [I] . 
Among them, the nitrogen-vacancy (NV) color center in diamond can be manipulated 
with full control and has remarkable properties for many applications ranging from 
high sensitivity and high resolution magnetometry [2jj6], to quantum information 
processing [7 



12 



and imaging in life science 13 15 . Most of these applications 



rely on the long coherence time of the NV defect electronic spin, which is mainly 
limited by magnetic interactions with a bath of paramagnetic impurities inside the 
diamond matrix. Therefore efforts have concentrated on creating diamond samples 



with a low concentration of impurities 16,17 , controlling the implantation of single NV 



defects 18 22 , manipulating and controlling the dynamics of their spin bath 23 27 



and understanding the interaction between the central spin and its environment both 



experimentally and theoretically 28 34 



The decoherence of a central spin in the presence of a spin bath has been addressed 
using several approaches [35| - [40] . Here we study both experimentally and theoretically 
the free induction decay (FID) of the electronic spin associated with a single NV defect in 
diamond. In particular, we have performed a statistical study on the coherence time, 
as a function of the strength of a magnetic field oriented along the NV defect symmetry 
axis. Our results indicate an increment of the coherence time at large magnetic fields 
in agreement with our numerical simulations and theoretical studies. In addition, we 
study the coherence time for different concentrations of the spin bath and identify the 
main features of the central spin dynamics. By dividing the bath into shells and cones, 
we analyse the contribution to decoherence of each impurity with respect to its position 
relative to the central spin. These results, which complement recent works reported 
in Refs 34 , 41 , might be particularly useful for diamond based magnetometry of dc 



fields (242 , diamond-based quantum information experiments such as nuclear spin bath 
control and nuclear-spin based quantum memories studies with isotopically engineered 



diamond samples 12 



The paper is organized as follows. We first introduce in Section 2.1 a theoretical 
model of a central spin immersed in a spin bath, which allows us to infer the phase 



memory time T 2 * of the NV defect electron spin. We then describe in Section 2.2| the 



numerical and experimental methods used to analyze the FID signals recorded from 
single NV defects. The experimental and numerical results are finally presented in 
Section |3j and compared to the theoretical predictions. 

2. Model and methods 



2.1. Theoretical model 

We consider the electronic spin S = 1 associated with a single negatively-charged NV 
defect in diamond (figure pla)). This central spin is nestled in a high-purity diamond 
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Figure 1. (a) Atomic structure of the NV defect in diamond consisting of a 
substitutional nitrogen atom (N) associated with a vacancy (V) in an adjacent lattice 
site of the diamond matrix. The NV defect axis z provides an intrinsic quantization axis 
for the electron spin and a magnetic field B is applied along this axis, (b) Schematic 
view of a single NV defect (red arrow) placed in a nuclear spin bath. Blue arrows 
indicate 13 C nuclear spins located randomly in the diamond lattice. Decoherence of 
the central electronic spin is induced by hyperfine coupling with the 13 C nuclear spins. 



lattice, where electron spin impurities like nitrogen donors (PI centers) are below 1 ppb 
and thus do not contribute to the decoherence of the central spin. Each lattice position 
can be occupied either by 12 C atoms (spinless) or by 13 C isotopes (nuclear spin I = |), 
which form a nuclear-spin bath (figure [ljb)). The natural abundance of 13 C in diamond 
is Pnat = 1.1% while isotopically-modified diamond samples exhibit 13 C concentrations 



as low as p = 0.01% 43 . In a high purity diamond sample, the decoherence of the NV 
defect electron spin is dominated by hyperfine coupling with the 13 C nuclear spins. The 
full Hamiltonian of the system is given by 

H = iJ NV + H hsLth + i^NV-bath • (1) 

In a magnetic field B, the Hamiltonian describing the NV defect H^y, the bath 
ifbath and their interaction i?NV-bath read 

H^ = DS 2 z + le B-S, (2) 
tfbath = E 7nB • I (n) + £ I (n) • ■ l(' m \ (3) 

n n<m 

#NV-bath = E S - AW - l(n) > ( 4 ) 

n 

where D/2n = 2.87 GHz is the zero field splitting of the NV defect electron spin, 
7 e /27r = 2.8 MHz/G (resp. 7 n /27r = 1.07 kHz/G) is the gyromagnetic ratio of the 
electronic (resp. nuclear) spin, is the hyperfine tensor between the central spin and 
nuclear spin n, and & nm ) is the dipole-dipole interaction between nuclei n and m. We 
note that we consider magnetic field amplitudes large enough to neglect strain-induced 
splitting terms in equation (2). The decoherence of the NV defect electron spin in the 
very low magnetic field regime (B < 0.1 G), for which strain- induced splitting protects 
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the central spin against magnetic field fluctuations, has been addressed elsewhere 44 . 

We focus on the free induction decay (FID) of a coherent superposition of the central 
spin as a function of the amplitude of a magnetic field applied along the NV defect axis, 
denoted z (figure [jja)). In this condition, the electron spin quantization axis is fixed 
by the NV defect axis and the eigenstates of H^y are denoted as \m s ) with m s = 0, ±1. 
The FID signal can be obtained by performing a Ramsey sequence | — r — | on the 
central spin, where a first | pulse rotates the central spin from a prepared state |0) to 
a coherent superposition |+) = ^ (|0) + |1)) which evolves over a time r. During such 
a free precession time, the central spin only interacts with the external magnetic field 
and the nuclear spin bath. Finally, a second | pulse rotates the central spin again, thus 
converting residual spin coherence into population. The FID signal can be interpreted 
as the probability to retrieve the initial state |0) and can be formally calculated as 

s(r) =Tr bat ,[(0|p(r)|0>], (5) 
p(r) = U T (T)p e g> pbathU T (Ty , (6) 

where the trace is taken over the Hilbert space of the bath. Here p e = |0)(0| is the 
initial electronic spin density matrix and Pbath — I\ n Pn ls the initial density matrix of 
the nuclear spin bath. The total unitary evolution operator of the system is denoted 
Ut{t) = -R(|)t/(r)i?(|), where -R(f) corresponds to a |-rotation operator over the 
electronic spin and U{t) = exp(— iHr). 

We neglect the slow dynamics related to dipolar interactions between nuclear spins 
within the bath, which is at most 2 kHz for nearest neighbor interaction. Indeed, 
experiments have shown that the NV defect coherence time T 2 * does not exceed few /is 
in diamond samples with natural abundance of 13 C 31 , indicating that the interaction 



between the central spin and the bath induces decoherence at a much faster rate than 
the dipolar coupling between 13 C's within the bath. In addition, the large value of the 
zero field splitting allows us to perform the secular approximation for the electronic 
spin even at small magnetic field values. Therefore, the terms proportional to S x , S y 
are neglected in the Hamiltonian. For a magnetic field pointing along the z-axis, the 
Hamiltonian can thus be written as 

H = DS 2 Z + le BS z + J2 7nS/i n) + E S Z A^ , (7) 

n n 

where the last two terms sum over each nuclear spin n and A^ n ^ = (A^ , A^j , A^) is 
now a hyperfine vector. For Ramsey experiments, we focus on the rotating frame at 
which the electron spin transition |0) — > |1) is addressed so that the Hamiltonian reduces 
to the last two terms in equation ([7]). We note that this Hamiltonian is diagonal in 
the electronic spin subspace and therefore it can be considered as a Hamiltonian for the 
nuclei, H mg , that is conditional to the state of the electron spin m s . The full Hamiltonian 
of the system can thus be written as 
i 

H= E \m s )(m s \ ® H ms , (8) 

m s =0 



Free induction decay of single spins in diamond 



5 



fl"m. = E fl S = E n S- l(B) . ( 9 ) 

n n 

where flj^ = 7„B + m s A' n ' represents the vector around which the nuclear spin n 
rotates, i.e. the Larmor vector. The evolution of the bath can then be written as the 
direct product of the evolution operator of each nuclear spin, U ms (r) = n n U^(r), 
where U$(t) = exp(-ztfgr). 

It follows from equation ^ that the FID signal is given by 

<r) = \ - ^Ren<Sn(r) , (10) 

where the product runs over all nuclear spins and S n = Tr„ p n UiUQ is the 
contribution of nuclear spin n to the signal. In order to obtain the effect of the nuclear 
spin bath on the FID signal envelope, we exclude the contribution of the intrinsic 14 N 
nuclear spin of the NV defect (1 = 1) and of strongly coupled nuclear spins because they 



only contribute to coherent oscillations of the signal 29 . For all other nuclear spins, we 



assume an initial state given by the high temperature limit, p n = 1/2. In this case, 

q(i) (n) (») (n) 

SJt) = cos — - — cos — h cos B {n> sin — - — sin — — , (11) 

where /3^ n ^ is the angle between fi^ and f2^ n ' > (see Appendix). 

As a figure of merit we consider the envelope of the FID signal, 

E(r) = ReY[S n (r) , (12) 

n 

where the product runs over lattice positions occupied by a 13 C atom. For a diamond 
lattice with a given concentration p of 13 C not all positions of the lattice are occupied 
by a 13 C. We therefore define the random variable x n = {1, 0}, with mean p, to denote 
the presence (x n = 1) or absence (x n = 0) of a 13 C at the lattice position n. For a given 
configuration of 13 C {x n }, the FID envelope can be written as 

E(t, {x n }) = Re]J[x n S n (r) + 1 - x n \ . (13) 

n 

We note that the product now runs over all positions of the lattice. Assuming that 
random variables independent, the average envelope is, 

E(r) = Rel[\pS n (T) + l-p] . (14) 

n 

We will use this last equation to analyze the average trend of the FID signal in the 
next sections. We now briefly discuss two simple limits to gain an insight into the FID 
dynamics. 

In the weak magnetic field limit (j n B A^), the FID envelope is given by, 

^(T)«nfl-> (B H a )=e- (T ^ )a , (15) 



f|E[^ (n) ] 2 l , (16) 
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where [A^] 2 = [A^] 2 + [A^>] 2 + [A£>] 2 . The decay of the FID signal is Gaussian and 
decoherence is caused because each nuclear spin precesses differently with magnitude 

AH 

In the large magnetic field limit, ^ n B ^> A^ n \ since cos (3 {n) ~ 1, the FID envelope 
reads 



E(r) 



1 2,HF 



n ■ , 

n x u 

(flM>]') 



,-(-/r 2 * HF ) 2 



(17) 



-1/2 



The signal also exhibits a Gaussian decay and decoherence is only due to the difference 
in Larmor precession along the direction of the external magnetic field because the 
anisotropic components of the hyperfine interaction are suppressed by the nuclear 
Zeeman energy. Since A^ 1 ) < A^ is always true, the NV defect coherence time is 
enhanced at high magnetic field. The enhancement factor 77 can be calculated as 



n 



J 2,HF 
2 2,LF 



(19) 



As justified in section 2.2 we neglect the contact interaction and therefore consider a 
pure dipole-dipole interaction between the central spin and the nuclei. The hyperfine 
components of nuclear spin n can then be written as 

At ] = a„,(l-3cos 2 ^ 



+ 



A 



3a n sin 9 n cos 9, n 



(20) 
(21) 



A (n) 



A 



(n) 



+ 



A(n) 



1 



3cos 2 ft„. (22) 

where a n is the magnitude of the interaction and 9 n is the angle between the NV 
symmetry axis and the imaginary line that joins the central spin and nucleus n. 
Considering a homogeneous distribution of nuclei and integrating over 9 and 0, the 
average enhancement factor is finally given by 

i7 = ^ = t/5 = 1.58..., 
T* LF V2 

independent of the concentration of nuclear spin impurities. 



(23) 



2.2. Numerical and experimental methods 

The FID signal is numerically simulated using two different approaches. The first 
method (Nl) consists in a statistical analysis over a large number of randomly sorted 



13 C distributions {x n } for which equation (13) is used to calculate the envelope of the 
FID signal. For each 13 C distribution, the result of the simulation is fitted to the 
function exp [— (t/T^ ) e ] in order to extract the exponent of the FID signal decay e and 
the coherence time T£. The second approach (N2) consists in calculating the average 



envelope of the FID signal using equation (14), thus considering that each lattice site has 
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a probability p to host a 13 C. This approach takes less computational time to infer the 
mean behavior of the single electronic central spin which otherwise must be obtained 
using the first approach by averaging the signal for a large number of different 13 C 
distributions. For both methods, a lattice of 100,000 atoms is considered, corresponding 
to approximatively 1000 nuclear spins for a diamond sample with natural abundance of 
13 C (p n at = 1-1%). In addition, the nuclear spin bath is divided into shells and cones 
in order to study the contribution to the FID dynamics of each 13 C with respect of its 
position relative to the central spin. 

We only consider the dipolar part of the hyperfine interaction in our model. 
Although a strong contact interaction has been measured for several 13 C's close to 
the central spin 45 , 46 and calculated theoretically 47 , we restrict the study to low 
13 C concentrations p < 4%, where decoherence is dominated by the dipolar interaction 
with the nuclear spin bath 31 . In this framework, strongly coupled nuclear spins only 
contribute to coherent oscillations of the FID signal but do not change its envelope. 
Therefore, we exclude from our simulations the lattice sites associated with nuclei that 
strongly interact with the central spin such as the well known 130-MHz-splitting linked 
to the nearest neighbor sites of the vacancy 48 and those with a contact interaction 
larger than few MHz. The neglected lattice positions correspond to about 100 lattice 
points, i.e. on average one nuclear spin for a natural abundance of 13 C. 

From the experimental side, individual NV defects hosted in a high-purity diamond 
crystal with a natural abundance of 13 C are optically isolated at room temperature using 
a scanning confocal microscope. A laser operating at 532 nm wavelength is focused 
onto the sample through a high numerical aperture oil-immersion microscope objective. 
The NV defect photoluminescence (PL) is collected by the same objective and directed 
to a photon counting detection system. Under optical illumination, the NV defect is 



efficiently polarized into the m s = spin sublevel 49 . This process provides an efficient 



state preparation for the Ramsey sequence. In addition, the PL intensity is significantly 
higher (~ 30%) when the m s = state is populated. Such a spin-dependent PL response 
enables the detection of electron spin resonances (ESR) on a single NV defect by optical 
means 50 . As discussed in the previous section, the FID signal is measured by applying 
a Ramsey sequence | — r — | to the NV defect electron spin (figure |2^a)). For that 
purpose, the ESR transition |0) -H- |1) is driven with a microwave field applied through 
a copper microwire directly spanned on the diamond surface. In addition, a permanent 
magnet placed on a three-axis translation stage is used to apply a static magnetic field 
with controlled orientation and amplitude. The magnetic field is aligned along the NV 
axis with a precision better than 1°. 
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3. Results and Discussion 

3.1. Enhancement of the coherence time at high magnetic field 

We first study the dynamics of the FID signal as a function of the amplitude of a 
magnetic field applied along the NV defect axis for a natural abundance of 13 C. Typical 
FID signals recorded from a single NV defect at different magnetic field magnitudes are 
shown in figure [2^b). Since the intrinsic nitrogen atom of the NV defect is a 14 N isotope 
(99.6% abundance), corresponding to a nuclear spin 1 = 1, each electron spin state 
is split into three sublevels by hyperfine interaction. The ESR spectrum thus exhibits 



three hyperfine lines, splitted by = —2.16 MHz 48 and corresponding to the three 
nuclear spin projections. As a result, beating frequencies can be observed in the FID 
signal at low field (figure W(b), upper and middle traces). By increasing the magnetic 



field, the NV defect gets closer to the excited-state level anti-crossing (LAC) which 



occurs at -Blac ~ 510 G 51,52 . In this case, electron-nuclear-spin flip-flops mediated 
by hyperfine interaction in the excited-state lead to an efficient polarization of the 14 N 
nuclear spin [53] . Consequently, the beating frequencies gradually disappear in the FID 
signal while increasing the magnetic field magnitude and the oscillations are then only 
due to the detuning of the microwave excitation (figure [2^b), lower trace). 

From these measurements, the exponent of the FID decay e and the 
coherence time T% are extracted through data fitting with the function 
exp[— (t/T 2 *) € ] J2m=-i Pm cos[27r(<5 + mA N )r], where 5 is the microwave detuning and 
(3 m is the population of each 14 N nuclear spin state. As expected from the model de- 
veloped in the previous section, we observe a significant enhancement of the coherence 
time with the amplitude of the magnetic field (figure [2^b)). In the high field limit, 
since the nuclear Zeeman energy is much larger than the strength of the hyperfine cou- 
pling, the quantization axis is fixed by the external magnetic field for all 13 C nuclear 
spins and the anisotropic components of the hyperfine interaction A ± are suppressed 



(see equation (18)). In the weak field limit, the quantization axis of each nuclear spin is 
rather given by their hyperfine fields and all the components of the hyperfine interaction 
contribute to decoherence of the central spin (see equation (Jl6|)). For a specific distribu- 
tion of 13 C, the numerical simulation reproduces well the behavior of the experimental 
results, as shown in figure W[c). 

To gain more insights into the FID dynamics, the coherence time T| and the 
exponent of the FID decay e are plotted in figures [3]^a) and (b) as a function of the 
magnetic field amplitude for two different NV defects. For both, the coherence time 
increases with the magnetic field strength. In addition, we also observe that the exponent 
of the FID decay deviates from e = 2 in the magnetic field range where the coherence 



time is rising, as recently reported in Ref. 41 . In this intermediate magnetic field 
regime, the nuclear Zeeman energy and the hyperfine interaction are of the same order 
of magnitude. Dephasing of the central spin is due to comparable nuclear rotations along 
non-parallel precession axis, fio and Q\. A more detailed expression for the contribution 
from a single nucleus to the signal is given on the Appendix for the intermediate regime. 
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Figure 2. (a) Experimental sequence used to measure the FID of the NV defect 
electron spin. A first laser pulse (2 /is) is used to polarize the electron spin in the 
m s — sublevel. A Ramsey sequence consisting of two resonant microwave (MW) tt/2- 
pulses separated by a variable free evolution duration r is then applied and a second 
laser pulse is finally used for spin state read-out. For data analysis, the NV defect PL 
recorded during the first 300 ns of the laser pulses is used for spin-state read-out while 
the PL recorded during the last 300 ns is used as reference, (b) Typical FID signals 
recorded from a single NV defect (NV12) for three different magnetic field amplitudes 
applied along the NV defect axis. The red solid line is data fitting as explained in the 
main text. The results of the fitting for the different magnetic fields are as follows. 
For B = 19 G, /3 +1 = 0.347, /3 = 0.333, = 0.320 and S = -0.61 MHz. For B = 77 
G, /3+i = 0.565, p Q = 0.288, /3_i = 0.147 and S = 3.65 MHz. And for B = 340 G, 
P+i = 1, /3 = 0, P-i = and S = 3.63 MHz. (c) Numerical simulations of the FID 
signal for a particular 13 C distribution. These graphs are obtained by multiplying the 
simulated FID envelope by the magnetic field dependent populations /3 m of each 14 N 
nuclear spin state. 



The simulation of the FID signal using the numerical method Nl agrees fairly well with 
our experimental measurements, as shown in figures^c) and (d). The average coherence 
time over 1000 different configurations of 13 C is shown on figure |3^c) as a function of 
the magnetic field. The coherence time increases around B = 150 G and at the same 
time the exponent of the exponential fit deviates from a Gaussian decay (e = 2) (Fig. [3] 
(d)). 

In order to extract a quantitative information of the enhancement factor rj, the 
coherence time was measured at low field (B m 20 G) and at high field (B rs 500 G) 
for a set of 20 single NV defects. Figure Hi a) shows a histogram of the coherence time 
in both regimes. A large variation of dephasing times is observed owing to the random 
distribution of 13 C around the central spin. The mean decoherence time at low field is 
(T 2 * (20 G)) cxp = 2.3 ± 1 fis while in the high field limit (T 2 * (500 G)) cxp = 3.7 ± 2 fjs 
(errors are standard deviations). The histogram of the enhancement factor 77 is shown 
in figure [41(b) , corresponding to (77) = 1.7 ± 0.5, in good quantitative agreement with 



the model developed in the previous section (see equation (23)). The histogram of the 
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Figure 3. (a) Coherence time T| and (b) exponent of the FID decay e as a function of 
the magnetic field strength for two single NV defects (NV12 and NV17). (c) Average 
coherence time and (d) exponent decay extracted from the numerical simulation of the 
FID signal over 1000 different nuclear spin bath configurations with p = 1.1%. Error 
bars indicate standard deviation. 
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Figure 4. (a) Experimental and (c) theoretical histograms of the coherence time T£ 
over 20 and 1000 different 13 C bath configurations, respectively, (b) Corresponding 
experimental and (d) theoretical distributions of the coherence time enhancement 
between weak magnetic field (B = 20 G) and high magnetic field (B = 500 G). 
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Figure 5. (a) Coherence time as a function of the magnetic field strength for four 
different concentrations of 13 C nuclear spin impurities, (b) Coherence time normalized 
to its value at high field (T^hf) as a function of the magnetic field magnitude for 
different concentrations of the nuclear spin bath. The enhanced factor is roughly 
v/5/2 = 1/0.632, independent of the concentration of impurities, (c) Coherence time 
as a function of the 13 C concentration p. The simulations agree with the dash-dotted 
line proportional to 1/y/p in the limit of either a weak or a high magnetic field. 



coherence time extracted from the numerical simulation over 1000 different nuclear spin 
bath configurations is shown in figure [4(c), leading to (T 2 * (20 G)) N i = 2.4 ± 0.7 /is and 
(T| (500 G))ni = 3.9 ± 1.4 /xs, meanwhile the histogram of the enhancement factor is 
shown in figure Bd), leading to (^)ni = I- 7 ± 0.4, in agreement with the experimental 
data. 

We note that the spread of coherence times increases at high magnetic field 
indicating that for some particular distributions of the nuclear spin bath there is a 
sizable enhancement, meanwhile for others there is little enhancement, as shown in 
figure [3(a) for two single NV defect having similar coherence times at low magnetic field. 
This feature is further illustrated in figures [4(c) and (d) which show the histograms of 
the enhancement factor inferred from the experiments and the numerical simulations, 
respectively. The contribution to the coherence time enhancement of each 13 C with 



respect to its position relative to the central spin is addressed in Section |3.3| where the 
bath is divided into shells and cones. 



3.2. Effect of the C nuclear spin concentration 

We now study the dynamics of the FID signal as a function of the magnetic field 
amplitude for different concentrations of 13 C nuclear spins p. For this purpose, we 
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use the numerical method N2 which directly provides the average trend of the FID 
signal within short computational time. Figure [5|a) shows the coherence time extracted 



from the envelope of the FID signal given by equation (14), as a function of the 
magnetic field strength for four different concentrations of impurities. By decreasing 
the content of 13 C, the coherence time of the central spin gets longer because 13 C atoms 
are placed on average farther away from the central spin, and thus exhibit weaker 



hyperfine interaction 17 . In addition, the magnetic field amplitude required to freeze 
the evolution due to the anisotropic component of the hyperfine interaction gets larger as 
the concentration of 13 C nuclear spins increases. As a result, the magnetic field strength 
linked to the transition between the coherence time in the low field limit (T£ LF ) and in 
the high field limit (T 2 * HF ) decreases with the 13 C content. We note that for a given 13 C 
content, the inner distribution of nearby nuclei might resemble to a bath configuration 
corresponding to smaller concentration. For this particular case, we expect to have a 
larger T 2 and a smaller magnetic field amplitude to increase T 2 , as it might be the case 
for NV12 in figure ||a). 

The coherence time normalized to its value in the high field limit is shown in 
figure [5]^b). The enhancement of the coherence time matches very well the value given 



by equation (23), T^hp/^lf = = 1/0.632, and is independent of the impurity 

concentration. Finally, the coherence time is plotted as a function of p on figure ^c). 
In the limit of either a weak or a high magnetic field, the coherence time scales as 1/ y/p, 



as expected from equations (16) and (18). In intermediate regimes, the scaling deviates 
from this simple behavior owing to the increment of the coherence time starting when 
7 n 5 ~ (see figure ga)). 



3. 3. Partition of the nuclear spin bath into shells and cones 

In this section, we analyse the contribution to decoherence of each nuclear spin impurity 
with respect to its position relative to the central spin, using the numerical method N2 
for a natural abundance of 13 C. 

First, we group 13 C according to their radial distance to the central spin. The 
nuclear spin bath is thus divided into shells Si with a 5 A width. The number of lattice 
sites in each shell is roughly proportional to the area of the shell times the width, as 
shown in figure [6]^a). To conduct the analysis, we infer the coherence time T 2S . from 
the FID envelope due to the nuclei belonging to the shell SV 

For low and high magnetic fields, the total coherence time can be calculated as 

2 = E(^M 2 - (24) 

i \ 1 2,Si J 




The coherence time for each shell at low and high magnetic fields is shown in figure |6j 
As expected, the closest shell to the central spin contributes the most to decoherence. We 
note that we have not considered strongly interacting nuclei, which results in neglecting 



the contribution of all lattice points of the first shell Si, as explained in Section J272 
We estimate the contribution to decoherence of each shell as (r 2 */T 2 * 5 J . Figure 6[c 
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Shell radius (A) Shell number Shell number 



Figure 6. (a) Shell partition of the nuclear spin bath. Each shell Si has a thickness of 
5 A. The number of lattice points per shell is indicated, (b) Coherence time versus shell 
number for low (10 G) and high (500 G) magnetic fields. The calculation is conducted 
with the numerical method N2 for a natural abundance of 13 C. The coherence time 
enhancement for each shell is equal to y/5/2, as shown in the inset, (c) Contribution 
of each shell to dccohcrcncc calculated as (TI/XJ s .) 2 . As expected, nuclei close to the 
central spin contribute the most. 



indicates that the second shell with w 650 lattice points, which corresponds to ~ 7 
nuclear spins at p nat = 1.1%, contributes in 85% to decoherence, meanwhile the third 
shell contributes only 10%. In addition, we note that, on average, the coherence time 
enhancement factor is not modified while changing the radial distance of the nuclear 
spins with respect to the central spin (see inset on figure [6|b)). 

We now analyze the contribution of each nuclear spin as a function of its angular 
orientation 6 with respect to the symmetry axis of the NV defect. For this purpose, the 
nuclear spin bath is divided into 20 angular cones, as shown in figure [7^ a). Each cone 
contains on average an equal number of lattice sites, so that their relative contributions 
to decoherence can be better compared. The coherence time T^ c . of the central spin is 
inferred by calculating the FID envelope due to the nuclei belonging to a given cone G{ 
of the nuclear spin bath. 

Figure [7](b) shows the FID envelope at weak and high magnetic fields for three 
different cones Ci, C5, and Cm, containing nuclei close to the north pole, to the magic 
angle 6 m ~ 54°, and to the equator, respectively. For cones C\ and C 10 , the FID signals 
are almost identical at weak and high magnetic fields, whereas for cone C5 there is a 
large enhancement of the coherence time at high magnetic field. The coherence time 
enhancement factor 77 is plotted as a function of the cone mid angle in figure ^c). This 
behavior can be understood as follows. Within our model, the hyperfine interaction 
is assumed to be purely dipolar between the central spin and each nuclear spin of 



the bath. In this framework, equations (20) and (21) indicate that the anisotropic 



component of the hyperfine vector A± is zero at the poles (6 = {0,7r}) and at the 
equator (6 = 7r/2), meanwhile the A zz component is zero at the polar magic angle 
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Figure 7. (a) Cone partition of the nuclear spin bath. Each cone contains 
approximately the same number of lattice sites (~ 5000). Cones C\, C5, and C10 
contain nuclear spin close to the north pole, to the magic angle 9 m ~ 54°, and to 
the equator, respectively, (b) Envelope of the FID signal for cones C\ (dashed lines), 
C5 (solid lines), and C10 (dashed-dot lines) at weak magnetic field (blue) and at high 
magnetic field (red) . (c) Coherence time enhancement factor 77 as a function of the cone 
mid angle. The maximum enhancement occurs for nuclei with a polar angle 9 close to 



the magic angle. The solid line is the theoretical prediction given by Equation (25 1 



(9m = cos _1 (l/v^3) ~ 54°). Therefore, since the anisotropic component is suppressed 
for large magnetic fields, the main contribution to the coherence time enhancement 
comes from the nuclei with axial component A zz m 0, i.e. placed at the magic angle 
9m ~ 54°. This effect is an analogous phenomenon to magic angle spinning in solid-state 



nuclear magnetic resonance spectroscopy 54 . Following equation (23 ), the enhancement 



factor as a function of the polar angle 9 can be written as 

VI + 3 cos 2 9 



n{9) 



T 2,Hf(0) 



(25) 



1 - 3 cos 2 9 

This formula accurately matches the enhancement factor inferred from the numerical 
simulations, as shown in figure |7j^c). 

From the partition of the bath into shells and cones, we conclude that nuclear spin 
baths with a large number of nuclei close to the magic angle and close to the central 
spin exhibit a large enhancement of their coherence time when the magnetic field is 
increased. This feature also explains the large variation of the enhancement factor 
shown in figures Qc) and (d), which correlates with the presence of nuclei close to the 
magic cone. 



4. Conclusions 



We have studied the decoherence of the electronic spin associated with a single NV 
defect in diamond placed in a 13 C nuclear spin bath. By recording the FID signal as a 
function of the strength of a magnetic field applied along the NV defect axis, we have 
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demonstrated an increment of the coherence time at high magnetic field, in agreement 
with theoretical studies. Numerical simulations of the FID signal indicate that the 
coherence time enhancement is independent of the 13 C concentration and is mainly due 
to nuclei positioned close to the vacancy and close to the polar magic angle (6 ~ 54°) 
with respect to the symmetry axis of the NV defect. This work points out interesting 
dynamics of a single spin placed in a spin bath, whose understanding is required for the 
development of diamond-based quantum information processing where single spins in 
diamond are used as robust solid-state qubits. 
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Appendix A. 

In this appendix, we show how to obtain the expressions for S n (r) and the envelope 
E(t) in the low and high magnetic field limits (equations ( 11 ),( 15 ) and (17)). We start 

(Al) 



from equation ( 10 ) and consider, 



where 



S n {r)=^ n {p n u{ n) H n) } 



U ma {r) = exp [-i(fl ma ■ I)r] = cos 



Here cr is the vector of Pauli spin matrices and u r 
has been omitted for simplicity. Consider 



Ulu 



cos — — COS 



-1{<T ■ Ui ) COS 



2 2 



t{(T ■ U ) COS 



- %\cr ■ u ms 

^m s /\^"m, 

- sin 



sm 



(A.2) 



The superscript n 



sm 



+ (cr • u )(cr • Ui) sm — — sm 



(A.3) 
(A.4) 



2 2 v u/v ' 2 2 

where (cr • uo) (cr • ui) = Uo • Ui + icr • Uo x Ui. For p n = 1/2, when the trace is taken 
over the Hilbert space of nucleus n, only those terms not proportional to a will survive 
because Tr„ (cr • u m .) = 0. Therefore, 



COS COS h 



S n (r) 

where we identify Uq™" 1 



u 



W „(n) 







• u 



sm sm 



(A.5) 



u 



(n) 



Note that in our construction tt 

as 

S n (r) = cos 1 - 1 > + (cos/3^ - 1) sin 



as the angle between the two Larmor vectors. 
'jnB, for all nuclei. Equation (A.5) can be arranged 



sm 



(A.6) 
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where the second term is zero for both the low and high field limits and it can be 
considered as the transition function between the two regimes. In the low field limit 
B « 0, 

2 



SJt) 



COS 



1 



1 / tt^T 



1 



(A.7) 



From equation (14), the envelope of the FID signal is, 

p(l-l[A^rfj+l-p 



E(t) « Ren n 



n„ i 



P [A (n) r] 2 



In the high field limit (>y n B > A (n) ), cos /3 (n) ~ 1, and 

,,)(«) n(") 

<S„.fr 



RJ COS ^— ^- PS 1 

2 8 



[4> 



Vl 2 



where we used = (A { 1 )2 + ( 7n S + A^) 2 ) 1/2 « 7ri 

7 n .B + Therefore, the envelope signal in this limit is 

1, 



Et) 



RelL 



n„(i 



(Ai 



(A.9) 

1/2 



(A.10) 
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